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ABSTRACT 


Propagation  of  large-amplitude  plane  waves  was  studied  in  rock 
specimens  from  1/2  to  12  inches  thick  using  in-material  piezoresistive 
stress  gages.  Sioux  quartzite,  charcoal  black  granite,  and  Coconino 
sandstone  were  examined.  Good  records  were  obtained  in  the  granite 
experiments,  allowing  separation  of  the  precursor  from  the  main  wave 
to  be  observed  at  five  successive  locations  in  a  single  experiment.  It 
is  concluded  that  time-dependent  effects  on  wave  propagation  in  quartzite 
and  granite  are  not  significant  over  these  distances.  Seven  experiments 
in  granite,  with  dry  and  water-filled  open  cracks,  displayed  the  effect 
of  the  crack  in  erasing  the  precursor.  The  results  show  that  the 


presence  of  a  crack  volume  of  a  percent  or  so  greatly  changes  the 
structure  of  the  precursor.  Only  a  few  records  were  obtained  in  the 
sandstone  experiments,  and  these  are  only  partially  understood. 
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FOREWORD 


This  report  expands  and  revises  our  previous  Final  Report,  DASA  2579. 
Since  that  report  was  issued,  we  have  conducted  five  more  experiments 
on  granite.  These  results  have  been  incorporated  into  the  previous 
report.  The  work  discussed  here  was  performed  between  May  1,  1969, and 
September  30,  1970, and  between  March  23,  1971,  and  October  15,  1971. 
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I  INTRODUCTION  AND  SUMMARY 

Calculations  of  stress  wave  propagation  and  seismic  coupling 
resulting  from  detonation  of  a  nuclear  or  chemical  explosive  in  an 
earth  material  require  a  realistic  model  of  the  response  of  the  earth 
material  to  an  extremely  wide  range  of  stress  loads,  ranging  from 
thousands  of  kilobars  near  the  explosive  source  to  fractions  of  a 
kilobar  at  great  distances  from  the  source. 

In  the  past  few  years,  laboratory  shock  measurements  have  been 
extensively  an  input  information  for  computer  calculations  of 
stress  wave  propagation  resulting  from  explosions  in  earth  materials 
Although  much  progress  has  been  made  toward  understanding  this  process, 
there  is  often  great  disagreement  between  field  measurements  and 
calculations  based  on  laboratory  measurements.  Part  of  the  explanation 
this  discrepancy  may  be  that  the  thin,  homogeneous  specimens  used  in 
laboratory  measurements  do  not  adequately  represent  field  conditions 
where  large-scale  fractures  and  joints  may  be  common.  Also,  if  cime- 
varying  effects  are  important,  the  laboratory  measurements  on  thin 
specimens  will  not  be  directly  applicable  to  field  studies  without  an 
adequate  model  on  which  to  base  long-time  extrapolation. 

In  the  Present  research  program,  experiments  were  conducted  to 
study  the  propagation  of  large-amplitude  plane  waves  in  larger  rock 
specimens  than  have  been  used  before.  The  aim  was  to  evaluate  the 
effects  of  open  fractures  and  the  importance  of  time-varying  effects 
Hard  rocks,  (Sioux  quartzite  and  charcoal  black  granite  from  Minnesota) 
and  a  porous  rock  (Coconino  sandstone,  both  dry  and  water-saturated) 
weie  studied.  Most  of  the  rock  specimens  were  several  inches  thick 
tke  largest  being  12  inches,  in  contrast  to  the  usual  laboratory  specimens, 
which  are  about  1/4-inch  thick.  A  list  of  shots  fired  in  the  program 
is  given  in  Table  I. 

The  first  six  experiments,  conducted  in  Sioux  quartzite,  were 
concerned  primarily  with  developing  suitable  techniques.  In-material 
manganin  stress  gages  were  used  to  record  stress-time  profiles  at  intervals 
Of  depth  in  the  rock.  The  large  specimens  presented  difficulties  in 
gage  emplacement  and  recording.  The  main  problem  was  piezoelectric  noise. 


1 


Table  I 


LIST  OF  SHOTS  FIRED 


Shot  No. 

Explosive  Train 

Specimen 

Comments 

1 

P-120  +  4"  TNT 

Sioux  quartzite 

Triggering  failed 

2 

P-120  +  4"  TNT 

Sioux  quartzite 

Noisy,  data  in 

Table  II 

3 

P-40 

Sioux  quartzite 

Technique  Shot, 

Table  II 

4 

P-40 

Sioux  quartzite 

Technique  Shot, 

Table  II 

5 

P-120  +  4"  TNT 

Sioux  quartzite 

Data,  Table  II 

6 

P-40 

Sioux  quartzite 

Technique  Shot, 

Table  II 

7 

P-80  +  2"  Baratol  + 
1/2"  A1  flyer 

Dry  Coconino 
sandstone 

Technique  Shot 

8 

P-80  +  2"  Baratol  + 
1/2"  A1  flyer 

Wet  Coconino 
sandstone 

Technique  Shot 

9 

P-80  +  2"  Baratol  + 
1/2"  A1  flyer 

Dry  Coconino 
sandstone 

Data,  Table  V 

10 

P-80  +  2"  Baratol  + 
1/2"  A1  flyer 

Wet  Coconino 
sandstone 

No  records 

11 

P-40 

Granite 

Technique  Shot 

12 

P-120  +  4"  Comp.  B 

Granite 

Data,  Table  m  and  IV 

13 

P-120  +  4"  Comp.  B. 

Granite 

Data,  Table  III 

14 

P-120  +  4"  Comp.  B 

Granite 

Data,  Table  uj 

15 

P-120  +  4"  Comp.  B 

Granite 

Data,  Table  m 

16 

Mousetrap  Plane  Wave 
Generator  +  2"  Comp.B 

Granite 

Data,  Tabic  m 

17 

P-120  +  4"  Comp.  B. 

v. 

Granite 

Data,  Table  J I I 

18 

P-120  +  4"  Comp.  B. 

Granite 

Data,  Table  III 

19 

P-120  +  4"  Comp.  B. 

Granite 

Data,  Table  III 

20 

P-120  +  4"  Comp.  B. 

Grani te 

Data,  Table  III 

21 

P-120  +  4"  Comp.  B. 

Granite 

Data,  Table  III 
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which  was  accentuated  by  the  long  leads  required  in  the  experiments. 

This  problem  was  largely  overcome  by  modifying  the  recording  system 
and  the  gage  design  to  minimize  common-mode  piezoelectric  pickup. 

A  total  of  six  experiments  were  conducted  in  Sioux  quartzite  with 
a  maximum  stress-wave  propagation-distance  of  4  inches.  When  an  attempt 
was  made  to  machine  larger  rock  slabs  from  a  block  of  quartzite,  we 
found  that  internal  cracks  were  spaced  close  enough  to  prevent  machining- 
even  1-foot-square  slabs.  These  cracks  were  not  at  all  apparent  from 
visual  inspection  of  the  outside  of  the  rough  block.  Therefore,  we 
substituted  charcoal  black  granite  supplied  by  the  Cold  Spring  Granite 
Company  in  Cold  Spring,  Minnesota. 

The  eleven  experiments  in  granite  yielded  the  most  complete  records. 
Except  for  one  preliminary  test  shot,  each  experiment  included  five 
manganin  stress  gages.  The  development  of  the  precursor  was  observed 
at  X-inch  intervals  to  a  propagation  distance  of  6  inches.  The  data 
indicate  that  time-dependent  phenomena  are  not  significant  at  this  scale. 
The  main  phenomenon  observed  is  the  separation  of  the  precursor  from  the 
main  wave  with  increasing  travel  distance.  The  differences  in  precursor 
amplitude  and  waveform  with  increasing  travel  distance  are  not  great  and 
are  attributed  to  experimental  scatter  caused  primarily  by  large  grain 
size  and  heterogeneous  composition  of  the  rock  at  the  scale  of  the  gage 
elemen  ts . 

Seven  experiments  with  wet  and  dry  open  cracks  in  granite  showed  vet’y 

interesting  results.  The  precursor,  which  is  separated  from  the  main 
wave  by  the  time  the  crack  is  encountered,  is,  in  effect,  pushed  back 
into  the  main  wave.  With  further  travel,  the  precursor  once  again  begins 
to  separate.  Computer  runs  using  SRI-PUFF  were  made  and  compared  with 

the  experimental  records.  The  latter  have  wider  variations  in  waveform 
and  time  separation  than  the  computed  profiles.  The  results  suggest  that 
in  the  field  a  crack  volume  of  a  percent  or  so  may  greatly  change  the 
precursor  structure  of  a  shock  wave. 
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'  - UBl»een  trim  taoout  U.UU3-inch)  layers  of  Epon.  The 

plate-impact  loading  used  was  expected  to  produce  a  flat-topped  pulse,  but 
the  waveforms  showed  stress  decay  immediately  behind  an  early  peak. 

Further  experiments  would  be  necessary ‘to  verify  that  a  material  property 
of  the  rock  was  being  observed. 

II  EXPERIMENTAL  RESULTS 


A.  Sioux  Quartzite 

Sioux  quartzite  was  initially  chose;,  to  be  the  main  rock  used  in  the 
program.  It  is  monomi neral ic  (~  99%  quartz)  with  grains  averaging  about 
0.1  mm  in  diameter  (Ahrens  and  Gregson,  1964).  For  our  specimens,  the 
density  average  2.624  ±0.003  g/cm3  p„r  eight  specimens,  the  ultrasonic 
velocity  was  5.30  ±0.05  hm/sec,  but  tor  the  ninth  specimen,  the  velocity 
averaged  4.33  km/sec.  Six  experiments  were  completed  before  it  was 
decided  that  the  rock  was  not  suitable  for  the  size  of  experiments 
contemplated.  Internal  cracks  limited  the  size  of  machined  specimens 
that  could  be  obtained. 

Three  experiments  involved  machined  slabs  of  dry  Sioux  quartzite 
explosively  loaded  with  P-120  plane-wave  lenses  and  4-inch-thick  pads  of 
TNT.  In  the  first  of  these  shots,  timing  from  ionization  gages  at  the 
explosive-rock  interface  did  not  work,  so  no  records  were  obtained.  In 
the  second  shot,  timing  was  by  electronic  delay  from  the  firing  pulse. 

The  stress  gages  used  on  the  first  two  shots  had  sensitive  elements 

0.625  inch  long  made  oi  flattened  manganin  wire  0.003-inch  thick.  The  four 

leads  on  each  gage  were  of  0.001-inch  copper  foil  cut  0.063-inch  wide. 

The  records  obtained  on  the  second  shot  had  so  much  noise  pickup  from  the 
piezoelectric  quartz  grains  along  the  6-  to  8-inch-long  foil  leads  that 
the  records  could  not  be  in terpreted • 

Two  technique  shots  were  designed  to  investigate  the  cause  of  the 
noise,  each  using  two  stress  gages  sandwiched  between  4  x  4  x  1/2-inch 
slabs  of  Sioux  quartzite.  Explosive  loading  was  by  a  P-40  plane-wave  lens. 


*  Epon (R)  adhesive  911F,  Shell  Chemical  Company. 
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XX  V, 


,  iaLLtnea  to  aoout  0.0008  inch  thick  by  0.004 


inch  wide,  instead  of  the  copper  foil  used  previously.  Even  with  these 
narrower  and  shorter  leads,  pickup  from  the  piezoelectric  effect  in  the 
quartzite  was  still  too  high  to  permit  precise  readings  from  the  records. 
However,  we  noted  that  pickup  started  at  the  time  the  shock  hit  the  rock- 
explosive  interface  about  1/2  inch  from  the  gages  and  did  not  grow  much 
as  the  shock  progressed  past  the  gages.  concluded  that  the  radiation 
from  grains  at  distances  of  1/2  inch  and  greater  was  more  important  than 

we  had  previously  suspected,  and  therefore  much  of  the  electromagnetic 
pickup  would  be  common  mode. 

°"  the  next  shot  (No.  4),  one  gage  used  leads  emerging  fro,  the  side 
as  in  the  first  shot,  „hile  the  other  gage  used  wire  leads  e.erging  through 
holes  drilled  in  the  quartzite  slab  away  from  the  explosive.  A  Tektronix 
Type  0  differential  preamp  was  used  on  the  loads  of  the  gage  emerging  fro. 
the  side  of  the  shot.  Pickup  was  reduced  by  about  a  factor  of  5,  thus 
giving  acceptable  data.  The  gage  with  leads  emerging  from  the  back  of 
the  shot  showed  little  pickup  until  the  wave  progressed  past  the  gage, 
at  which  point  the  record  became  unacceptably  noisy,  beads  emerging 
from  the  back  of  the  shot  through  drilled  holes  are  both  more  expensive 
to  install  and  more  likely  to  interfere  with  the  wave  propagation  process, 
he  will  use  them  only  when  we  need  a  long  time  response,  since  they 

generally  do  not  break  as  quickly  as  gages  contained  completely  in  a 
crack  between  slabs. 

Shot  5,  on  12  x  12-inch  slabs  of  Sioux  quartzite  (Fig,  l)  used 
gages  with  1/4-inch-long  sensitive  elements  made  of  flattened 
0.002-inch  manganin  wire.  All  leads  were  of  flattened  0.002-inch 

copper' wire  emerging  from  the  corners  of  the  cracks.  The  voltage  leads 
were  laid  as  close  together  as  possible  (about  1  mm)  to  further  improve 
common  mode  rejection  of  piezoelectric  pickup.  The  stress  wave  record 
irom  gage  No.  1  is  shown  in  Fig.  2. 

The  transit  times  between  the  gages  at  2  inches  and  4  inches  were 
9.3  psec  for  the  precursors  and  9.9  ^sec  for  the  main  wave.  Average 
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PLANE  WAVE 
LENS 
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SIOUX 

QUARTZITE 


SURFACE  A 
SURFACE  B 


GAGE 
NO.  1 


GAGE  NO.  2 


SURFACE  A 
Looking  Toward  TNT 


GAGE  NO.  3 


SURFACE  B 
Looking  Toward  TNT 


GAGE  NO.  4 


GA-7942-1 

FIGURE  1  GAGE  ARRANGEMENT  FOR  SHOT  NO.  5 
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slab  before  the  shot  but  was  over  5.22  km/sec  on  all  other  slabs. 
Unfortunately  precursor  transit  times  were  not  measurable  for  the  other 
slabs.  An  unusually  low  ultrasonic  velocity  in  one  rock  specimen  is 
often  caused  by  an  unusually  high  number  of  microcracks.  A  high-amplitude 
precursor  will  close  the  cracks  so  that  precursor  velocities  have  less 
scatter  than  small  amplitude  ultrasonic  waves. 

Shot  6  was  done  to  check  whether  the  high-frequency  noise  on  the 
records  is  produced  by  stress  waves  from  formation  and  growth  of  cracks 
or  by  the  piezoelectric  effect.  Two  similar  stress  gages  were  built  into 
the  shot,  one  connected  to  our  usual  power  supply  and  the  other  with  no 
power  supply  attached.  Both  gages  should  have  picked  up  piezoelectric 
noise  equally,  but  only  the  gage  with  the  current  pulse  would  pick  up 
noise  from  fracture-produced  stress  waves.  Both  gages  picked  up  noise 
with  equal  RMS  amplitude;  hence  it  is  concluded  that  the  noise  is 
piezoel ectric . 

The  results  of  the  Sioux  quartzite  experiments  are  given  in  Table  II. 


Table  II 


SIOUX  QUARTZITE  EXPERIMENTS 


Note: 


Re«ults  In  parentheses  ire  regarded  as  poor  because  the 
Theae  data  .ere  omitted  from  calculations  of  means  and 


records  were  unusual  1 
standard  deviations. 


y 


noisy. 
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PRE¬ 

MAIN 

CURSOR  RAMp 

WAVE 

RISE  DURATION 

RISE 

TIME 

TIME 

GA-7942-3 

FIGURE  3  TERMINOLOGY  USED  TO  DESCRIBE  STRESS  WAVE 


Stress  amplitude  .as  calculated  Iron,  the  gaee  records  by  the  linear 
relationship  (Keough,  1968), 


_  AR  10 

°  13  Stl'eSS  1,1  kilobars»  ^  is  the  ratio  of  the  resistance 

change  to  the  initial  resistance  of  the  manganin  gage  element.  Since 

a  constant  current  source  was  used  and  the  resistance  and  resistance 
change  of  the  active  element  is  large  compared  with  that  of  the  leads  in 
the  speciman,  -  is  taken  to  he  equal  to  the  ratio  of  change  in  voltage 
to  initial  voltage  as  recorded  by  the  oscilloscopes. 

Shots  2  and  5  were  loaded  .1th  a  P-120  plane-wave  lens  plus  a  l-i„ch. 
thick,  12-i „eh-diameter  pad  or  TNT.  The  other-three  shots  were  loaded 
with  P-10  plane-wave  lenses  only.  Precursor  rise  time  and  amplitude, 
amplitude  at  the  back  „t  the  ramp,  and  main  .ave  rise  time  showed  no 
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«  a  uuuKness . 

'  °0n  thS  °"lv  nonsteac*y-state  phenomenon  observed  «es  the 

increase  an  r..p  length  „th  rock  thlckness  oauSed  bj.  tfce  ^ 

the  precursor.  Tab.e  „  gives  the  means  and  standard  deviations 
rise  times  of  the  precursor  and  main  »ave  and  tor  amp.itudes  tor  the 
precursor  at  the  trout  ot  the  ramp  and  at  the  bach  ot  the  ramp.  »ote 

**”  ““  "*  S1’llar  f°'  and  main  „ave.  A  measure 

f  °f  da‘a  IS  «*"■  *  “vidi„g  the  standard  deviation  by  the 

r,  °btaln  **  This  coetticient  is  smaller 

main  wave  rise  time  than  for  that  of  the  precursor.  The  coefficient 

J  ™°"  "  8130  “  ~  «  tf,e  of  the  ramp  than 

k,T  of  the  precursor  at  the  froM  ^  the  ^  conciude 

that  details  of  the  front  edge  of  the  precursor  are  more  sensitive  to 
roch  inhomogeneities  than  details  connected  with  the  rise  of  the  plane- 
wave  from  the  back  of  the  precursor. 

B.  Granite 

or  c  rru  Maok  era,,ite  suppiied  by  *he  c°id  spri,,e  c°»p*"y 

pring,  Minnesota,  was  substituted  for  Sioux  quartzite  when  the 
after  proved  to  be  unsuitable  for  machining  into  large,  uniform  slabs, 
our  granite  was  quarried  from  the  same  site  as  that  currently  being  used 
ngh  IITRI  in  a  project  involving  triaxial  testing  of 
very  large  samples.  The  rock  is  massive  with  no  noticeable  preferred 

orientation  of  grains  or  tpv(-i,nD  ,,  . 

•  Most  grains  are  between  1  and  5  mm 

“  diameter.  A  few  larger,  denser,  fine-grained  inclusions  are  present. 

These  showed  up  well  in  X  rays  of  1-  and  2-inch-thick  slabs,  allowing 

leceTfl°th0nly  "°re  U',lf0r"  SlabS  f°r  the  — auction  of  the  shots. 

Very  1<lrge  sPe=i»ans.  »hich  could  not  be  handled  easily 
an  specimens  were  routinely  checked  for  density,  ultrasonic  compressil.al 
velocity,  and  X-ray  uniformity.  Scatter  within 

.  .  scatter  within  a  specimen  and  from  specimen 

to  specimen  was  relatively  small.  Average  density  measured  was  2.719  * 

•  03  g/cm  ,  and  average  compressions!  velocity  measured  was  5.41  ±  0  12 
km/ sec.  x  ^ 

A  total  of  eleven  experiments  were  conducted  In  the  charcoal  black 

granite.  One  was  a  simple  test  shot  to  evaluate  piezoelectric  noise  and 
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experiments,  each  with  five  manganin  stress  gages.  In  Shots  12  and  13 
the  stress-time  profile  was  recorded  at  1-inch  intervals  in  6-inch-thick 
specimens.  In  Shots  14,  18,  20,  and  21  the  effect  of  a  water-filled 

inclined  crack  was  recorded.  In  Shots  15,  17,  and  19  the  crack  was  open 
and  dry.  In  Shot  16  gages  were  placed  at  2-inch  intervals  in  a  12-inch 
thick  specimen.  A  summary  of  results  is  given  in  Table  III. 


Shots  12  and  13  were  very  similar  except  for  the  manner  in  which 
the  gages  were  inserted  into  the  rock  assembly  (Fig.  4).  in  Shot  12, 
the  rock  was  cut  and  machined  into  slabs  1-inch  thick  and  more  than  12 
inches  across  all  diameters.  Six  of  these  slabs  were  stacked  to  assemble 
the  shot,  with  a  manganin  wire  gage  at  each  interlace.  The  leads  were 
brought  out  the  largest  dimension  to  the  sides.  The  shot  v.is  fired  with 
a  P-120  plane-wave  lens  and  a  4-inch  by  12-inch  diameter  cylinder  of 
Composition  B  high  explosive. 

In  shot  13  a  quarter-round,  6-inches  thick  and  more  than  12  inches 
across  all  diameters,  was  machined  flat  on  one  end  where  the  explosive 
was  to  be  placed.  A  2-1/2-inch  diameter  core  was  drilled  out  of  the 
center,  perpendicular  to  the  machined  surface.  One-inch-thick  discs  of 
granite  were  then  machined  to  fit  snugly  (<  0.010-inch  glue  line  across 
the  diameter)  back  into  the  hole.  Five  manganin  wire  stress  gages  were 
fitted  into  the  core,  one  at  the  center  of  each  disc  interface.  Copper 
leads  were  used  from  the  manganin  element  to  the  circumference  of  the 
core.  Leads  were  brought  out  to  the  back  of  the  shot  in  narrow  slots 
cut  into  the  circumference  of  the  core  and  filled  with  dental  amalgam. 

Good  quality  records  were  obtained  from  all  gages  for  shots  12  and 
13  (Table  III  and  Fig.  5).  The  two  shots  were  very  similar  in  showing  the 
rate  of  precursor  development.  The  precursor  is  slower  to  separate 
clearly  from  the  main  wave  in  granite  than  it  is  in  quart zl  I  o  ff'ig.  6) 
probably  because  of  the  larger  grain  size  and  the  more  varied  mineralogy, 
flic  scatter  in  precursor  amplitude  appears  to  be  random,  caused  by  large 
grain  size  and  heterogeneous  composition  of  the  rock  at  the  scale  of  the 
gage  elements.  The  precursor  amplitude  and  waveform  does  not  seem  to  be 
aflected  significantly  by  Lime  or  distance  of  propagation.  This  is  an 
impoitant  result  because  of  the  question  of  the  validity  of  small-scale 
laboratory  experiments  to  determine  rock  properties  that  are  required  for 
in  situ  problems. 
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Gage  Thickness 


26. 

,34 

78. 

,31 

78. 

,31 

76. 

,31 

104. 

37 

51 

.99 

76, 

.05 

78, 

.05 

104, 

.72 

104. 

CM 

51. 

,66 

78. 

,12 

78. 

,12 

103. 

,24 

103. 

24 

26. 

,95 

29. 

,40 

79. 

,40 

106. 

,65 

106. 

65 

Precursor 


Rise  Time 
(usoc^ 


Rise  Tine 

Ampl .  ot  Back 

(wee) 

1  (kbar) 

Main  Wove 


Ampl . 
(kbnr) 


1 

49.65 

1  Gage  failed  before  HEL  was  reached. 

2 

100,96 

0.4  1  26  1  1.3 

1  42 

.. 

3 

152.34 

0.6;  1  20  I  2.15 

40 

0.62 

140 

4 

J03.70 

No  records  obtalred 

5 

254.78 

No  records  obtained 

SHOT  17 

1 

2 

3 

4 

2 

26.44 

76.89 

78.89 

105. 

105.  •' 

0.82 

0.63 

0.44 

0.65 

35 

42 

22 

28 

0.41 

0.22 

0.30 

0.91 

0.93 

63 

42 

55 

60 

52 

0.14 

0.21 

0.24 

0.19 

0.26 

200 

178 

175 

172 

156 

SHOT  21 

31 

0.39 

165 

51 

0.41 

160 

56 

0.50 

170 

50 

0.28 

157 

42 

0.4G 

150 
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(b)  SHOT  NO.  13 
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EXPLOSIVE 
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Id)  SHOTS  NO.  17,  18,  AND  21 


PLANE  WAVE 
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EXPLOSIVE 
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le)  SHOTS  NO.  19  AND  20 


FIGURE  4  GAGE  ARRANGEMENTS  FOR  GRANITE  EXPERIMENTS  x3 
location  and  number. 
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indicates  gage 
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QUARTZITE 


FIGURE  6 


PRECURSOR  IN  GRANITE  COMPARED  WITH  PRECURSOR  IN  QUARTZITE 
BOTH  AT  2  INCH  DEPTH.  Shots  12  and  5. 


Gage  records  in  the  granite  experiments  were  good  enough,  in 
contrast  to  the  more  noisy  quartzite  records,  to  permit  the  use  of 
Lagrangian  analysis  to  derive  the  stress  particle  velocity  and  stress- 
volume  loading  curves.  Figures  7  and  8  and  Table  IV  show  average 
loading  curves  obtained  for  the  granite.  The  significance  of  time- 
dependent  effects  was  checked  by  doing  the  analysis  using,  first,  gages 
1,  2,  and  3,  and  then  using  gages  3,  4,  and  5.  The  resulting  curves  were 
nearly  identical,  as  shown  in  Fig.  8. 


FIGURE  7  STRESS-PARTICLE  VELOCITY  LOADING  CURVE  FOR  GRANITE. 
Obtained  by  Lagrangian  analysis  of  Shot  No.  12. 
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0.28  0.30  0.32  0.34 

VOLUME  —  cm3/g 

FIGURE  8  STRESS-VOLUME  LOADING  CURVES  FOR  GRANITE.  Obtained  by 
Lagrangian  analysis  of  Shot  No.  12. 


0.36  0.38 

GA-7942-12 


Lagrangian  analysis  is  described  in  general  by  Fowles  and 
Williams  (1970).  In  our  present  work  we  consider  that  stress  is  a 
function  of  particle  velocity  only,  which  simplifies  the  calculations. 
The  fact  that  the  loading  curves  derived  from  gages  1,  2,  and  3  are 
nearly  identical  to  those  derived  from  gages  3,  4,  and  5  indicates  that 
this  assumption  is  warranted.  The  computer  program  is  given  in  detail 
by  Murri  and  Smith  (1970). 
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Table  IV 


LOADING  CURVE  FOR  GRANITE, 

SHOT  NO.  12, 

GAGES  1,  3,  5 

Stress 

Particle  Velocity 

Volume 

Local  Sound 

(Jibar  ) 

(mm/  fJsec  ) 

(cm3/ g ) 

(mm/ /Usee  ) 

0 

0 

0. 3676 

5.905 

5 

0.0312 

0.3656 

5.858 

10 

0.0624 

0.3637 

5.806 

15 

0.0938 

0.3617 

5.759 

20 

0.1253 

0.3597 

5.709 

25 

0.1568 

0.3578 

5.661 

30 

0.1884 

0.3558 

5.618 

35 

0.2202 

0.3537 

5.545 

40 

0.2523 

0. 3517 

5.461 

45 

0.2849 

0.3495 

5.306 

50 

0.3180 

0.3474 

5.227 

55 

0.3513 . 

0.3451 

5.157 

60 

0.3848 

0.3429 

5.110 

65 

0.4184 

0.3406 

5.065 

70 

0.4521 

0.3384 

5.022 

75 

0.4858 

0. 3361 

4.987 

80 

0.5195 

0.3338 

4.950 

85 

0.5532 

0.3315 

4.915 

90 

0.5869 

0.3293 

4.882 

95 

0.6206 

0.3270 

4.850 

100 

0.6544 

0.3247 

4.818 

105 

0.6880 

0.3224 

4.785 

110 

0.7217 

0.3202 

4.753 

115 

0.7554 

0.3179 

4.721 

120 

0.7891 

0.3156 

4.689 

125 

0.8227 

0.3134 

4.657 

130 

0.8564 

0.3111 

4.622 

135 

0.8900 

0.3088 

4.587 

140 

0.9237 

0.3066 

4.550 

145 

0.9574 

0.3043 

4.513 

148 

0.9776 

0.3030 

4.491 
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In  the  calculations,  particle  velocity  is  first  obtained  by 


'i-1 


■kf 


dC 

dh/dt 


with  dh/dt  estimated  from  the  gage  data.  Here  u.  is  the  particle 


.  th 


velocity  at  the  i  point  (initially  u^  is  zero  for  compression 
calculations),  pQ  is  the  initial  density,  h  is  the  Lagrangian 
distance  coordinate  specifying  the  positions  of  the  gages,  or  is  the  stress, 
and  t  is  the  time.  Next,  the  volume  is  computed  by  integrating  the  equation 


V. 


=  V. 


where  V 


i-1 

is  the  volume  of  the 


-  —  (— 

0  2  J  (dh/ 


(dh/dt)' 


.  th 


point  (initially  V  is  the  zero 


s  tress  vol-umer 


i-1 


o 


The  local  sound  speed  is  then  calculated  from  the  equation 


V  Q 
i  dt 


where  C/  is  the  sound  speed  of  the  ith  point. 

Shots  14  and  15  were  designed  to  record  the  effects  of  open  cracks, 
dry  and  water-filled,  on  the  stress  profile.  In  both  experiments  a 
30°  angle  cut  was  made  in  a  2-inch-thick  slab  of  granite.  The  cut  was 
machined  and  held  open  0.005  inch  In  Shot  14  the  crack  was  filled  with 
water  under  vacuum  and  sealed;  in  Shot  15  the  crack  was  left  dry.  The 
assembly  is  shown  in  Fig.  4c.  The  explosive  chain,  a  P-120  and  a  4-inch 
pad  of  Composition  B,  is  in  contact  with  a  1-ir.ch  slab  of  granite,  which 
is  followed  in  turn  by  the  2-inch  slab  with  the  30°  cut,  a  l-inch  slab, 
and  a  covering  1-1/2-inch  slab.  Stress  gages  are  located  between  slabs 
of  granite.  Beneath  the  2-inch  slab,  two  gages  are  positioned  under 
the  open  cut  and  one  gage  on  the  other  side  of  the  cut. 

Cage  records  of  the  two  experiments  are  shown  in  Figs.  9  and  10. 

In  both  experiments,  the  records  from  the  gages  beneath  the  cuts  show 
an  erasing  of  the  precursor  as  can  be  seen  by  comparing  records  from 
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F'GURE  9  REC°RDS  FR°M  SH0T  G89eS  3  “«■  4  ™  belo*  water-filled 

crack.  kage  2  is  at  same  plane  but  not  below  crack. 


gnge  2  with  those  from  gages  3  and  4.  The  erasing  is  more  complete 
in  the  case  of  the  dry  crack  (Fig.  10). 

Computer  runs  using  SRI  PUFF  (Seaman,  1970)  have  been  -ade,  treating 
the  crack  as  a  layer  normal  to  the  direction  of  stress  wave  propagation. 

An  approximation  to  the  experimentally  measured  stress-time  pulse  was 
allowed  to  propagate  across  the  crack  to  the  recording  position  below 
the  crack.  The  three  cases  treated,  no  crack,  water-filled  crack,  and 
dry  crack,  are  shown  in  Fig.  11  along  with  experimentally  recorded  profiles. 
The  curve  from  Shot  15  in  Fig.  11b  has  been  adjusted  in  arrival  time  to 
agree  with  the  curves  from  Shot  14.  Several  important  observations 
can  be  made:  (1)  The  experimental  profiles  have  more  variations  in 
waveform  and  time  separation  than  the  computed  profiles.  This  may  be 
partly  experimental  scatter,  but  study  of  Figs.  9  and  10  shows  that 
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FIGURE  10 


GA  7942-9 


STRESS-TIME  RECORDS  FROM  SHOT  15.  Gages  3  and  4  are 
Gage  2  is  at  same  plane  but  not  below  crack. 


below  dry  crack. 


gage  records  3  and  4  from  below  the  crack  are,  in  each  case,  very  similar. 

n^f  Jh\rater"filled  CraCk  Pr°file’  as  co^Puted,  lies  closer  to  the  solid 
profile  than  to  the  dry  crack  profile.  The  experimental  records  show  the 

reverse.  However,  later  experiments,  discussed  below,  are  more  in  agreement 
with  the  computation.  (3)  The  computed  slopes  are  not  as  steep  as  the 
recorded  slopes.  This  is  a  fault  of  the  artificial  viscosity  computation 
and  could  be  improved  with  longer  runs. 

Shot  16  was  the  largest  experiment  of  the  program.  It  was  similar 
in  design  to  Shot  13  in  that  gages  were  placed  between  discs  of  rock  in 
a  cylindrical  core  that  was  inserted  into  a  machined  hole  in  the  center 
of  the  large  specimen.  Five  gages  were  used  at  intervals  of  2  inches. 

A  mouse-trap  plane-wave  generator  (Duvall  and  Fowles,  1963)  with  a 
geodesic  line  wave  initiator  (Erkman,  1959)  was  used  to  detonate  a 
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FIGURE  11  COMPUTED  AND  EXPERIMENTAL  GAGE  RECORDS  COMPARING  EFFECTS 
OF  WATER-FILLED  AND  DRY  CRACKS  IN  GRANITE 


24  x  24  x  2  inch  slab  oi  Composition  B  high  explosive  in  contact  with 

the  rock  surface.  An  experiment  was  conducted  in  an  attempt  to  sensitize 
nitromethane  so  that  it  could  be  used  as  the  explosive  to  be  initiated  by 

the  mouse-trap  arrangement  in  our  intended  experiments  on  2  x  2  ft  rock 
slabs.  With  a  mouse-trap  arrangement  we  attempted  to  initiate  plain 
nitromethane,  nitromethane  with  1%  triethylamine  added,  plain 
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GA-7942-10 


FIGURE  12  PRECURSOR  IN  GRANITE  AT  6  INCHES  FROM 

EXPLOSIVE-GRANITE  INTERFACE.  Shot  No.  16. 


nitromethane  with  floating  micro  balloons,  and  nitromothano  with  1", 
triethylamine  plus  floating  micro  balloons.  None  of  the  arrangements 
gave  a  satisfactory  detonation  from  the  mouse-trap  loading.  Therefore, 

»e  decided  to  use  the  more  expensive  and  less  uniform  machined  Composition 
B  slabs  for  the  explosive  in  the  large  shots. 


Results 
early  in  the 
recorded  the 
main  shock. 


for  Shot  16  are  summarised  in  Table  III.  Gage  1  shorted  out 
precursor  rise,  giving  only  time-of-arrival  data.  Gage  2 
precursor  development  and  failed  during  the  arrival  of  the 
Gage  3  shows  a  complete  record  of  the  shock  compression 
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(Fig.  12).  Unfortunately,  a  time  delay  generator  malfunctioned  causing 
the  loss  of  records  from  gages  4  and  5. 

Shots  17  and  18  were  identical  in  construction  except  that  the 
horizontal  0.006-inch  open  crack  was  dry  in  Shot  17  and  filled  with  water 
in  Shot  18.  Gage  records  are  shown  in  Fig.  13  and  14.  The  records  for 
gages  2  and  3, located  1  inch  beyond  the  crack,  show  that  the  dry  crack 
in  Shot  17  causes  the  precursor  to  have  a  steep  rise  to  about  30  kbar, 
whereas  in  Shot  18  the  precursor  rises  more  slowly  and  evenly  to  about 
60  kbar.  In  Shot  17  the  rise  from  0  to  50  kbar  takes  about  0.7  |isec;  in 
Shot  18  it  takes  1  M-sec.  In  the  first  case  the  crack  has  to  be  closed 
before  energy  is  transmitted  beyond  the  crack.  In  the  second  some  stress 
wave  energy  is  transmitted  directly  through  the  water.  A  series  of 
reverberations  permits  the  rest  of  the  energy  to  be  transmitted.  The 
result  is  seen  as  a  slower  rise  time  from  zero  stress  at  gages  2  and  3. 
Little  difference  is  seen  by  the  time  the  stress  wave  reaches  gages  4  and  5. 
Local  inhomogeneities  obscure  the  finer  details  of  the  waveform. 

o 

In  Shots  19  and  20  the  0.006-in;h  open  crack  was  inclined  at  60  to 
the  horizontal.  The  cracks  were  dry  in  Shot  19  and  filled  with  water  in 
Shot  20.  The  gage  records  for  the  two  shots,  shown  in  Figs.  15  and  16, 
are  not  significantly  different.  Presumably,  the  inclination  of  the  crack 
causes  the  effect  on  the  precursor  to  be  spread  out  over  time  and  over  the 
area  of  the  gage,  making  the  differences  between  the  wet  and  dry  cracks 
less  obvious. 

Shot  21  was  identical  to  Shot  18,  with  a  water-filled  horizontal 
crack,  except  that  in  Shot  21  the  crack  width  was  increased  to  0.024  inch 
from  the  0.006  inch  of  No.  18.  The  records  are  shown  in  Fig,  17.  In 
general  they  are  similar  to  those  of  Shot  18  (Fig.  14) .  However,  a 
pronounced  effect  on  the  shape  of  the  precursor  can  be  seen  in  Fig.  18, 
which  compares  gage  records  at  similar  positions  in  Shots  12  and  13,  17, 

18,  and  21.  In  each  case  the  average  of  two  gage  records  is  plotted  with 
times  adjusted  for  simultaneous  arrival  of  the  main  wave.  The  average 
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FIGURE  14  STRESS-TIME  RECORDS  FROM  SHOT  18 


FIGURE  15  STRESS-TIME  RECORDS  FROM  SHOT  19 


FIGURE  16  STRESS-TIME  RECORDS  FROM  SHOT  20 
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FIGURE  18  EXPERIMENTAL  RECORDS  OF  EFFECTS  OF  CRACKS  ON  STRUCTURE 
OF  PRECURSOR  IN  GRANITE 
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of  gages  2  and  3  for  Shot  18  shows  much  less  precursor  structure  but 
just  as  much  overall  risetime.  The  average  of  gages  2  and  3  in  Shot  21 
shows  almost  no  precursor  structure  and  less  rise  time.  The  similar 
record  for  Shot  17  with  the  dry  crack  is  shown  for  comparison.  As  noted 
above,  the  precursor  structure  remains,  but  the  risetime  is  reduced 
significantly.  Since  gages  2  and  3  are  1  inch  from  the  crack,  the  crack 
itself,  even  in  Shot  21,  is  only  2.4%  of  the  distance  from  the  crack  to  the 
gages.  In  Shot  18  the  crack  is  only  0.6%  of  that  distance.  These  results 
suggest  that  in  the  field  a  crack  volume  of  a  percent  or  so  may  greatly 
alter  the  precursor  structure  of  a  shock  wave. 

C .  Coconino  Sandstone 

Coconino  sandstone  was  chosen  for  use  as  the  porous  rock  in  the 
program.  It  is  nearly  monomineralic  (97%  quartz,  3%  feldspar,  with 
a  weak  silica  cement).  The  grain  size  is  0.1  to  0.2  mm  in  diameter 
(Ahrens  and  Gregson,  1964).  Our  specimens  had  a  porosity  of  about  20%. 

The  ultrasonic  sound  velocity  was  3.27  km/sec. 

Four  experiments  were  performed  on  Coconino  sandstone,  two  with  the 
rock  dry  and  two  with  the  rock  saturated  with  water.  The  first  two 
shots,  one  wet  and  one  dry,  were  otherwise  similar,  consisting  of  two 
slabs  of  sandstone,  approximately  8  x  10  x  l/2-inch  thick,  cut  parallel 
to  the  bedding.  Three  gages  were  inserted  in  the  crack  between  the  two 
sandstone  slabs,  each  gage  in  a  different  insulating  material.  The  main 
purpose  of  these  preliminary  experiments  was  to  check  the  effectiveness 
of  glass,  mica,  and  Epon  adhesive  as  insulating  and  armoring  material 
to  prevent  shorting  or  breaking  of  in-material  gages  in  the  porous 
sandstone.  An  additional  manganin  gage  was  used  on  the  back  surface  of 
each  shot  in  a  C-7  buffer. 


*  Epon  (R) ,  adhesive  911F,  Shell  Chemical  Company, 
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The  loading  of  each  shot  was  by  impact  of  a  1/2-inch  aluminum  flyer 
Plate  driven  by  a  P-80  wave  lens  plus  a  2-inch  Baratol  pad.  The  glass- 
insulated  gages  did  not  work  in  either  shot;  traces  were  too  noisy  to 

read.  Mica  and  Epon  insulation  worked  in  both  shots,  although  the  Epon 

insulation  gave  a  longer  recording  time.  Mica-insulated  gages  indicated 
114-kbar  peak  stress  in  the  dry  shot  and  160-kbar  in  the  wet  shot.  Epon- 
insulated  gages  indicated  128-kbar  peak  stress  in  the  dry  shot  and  174- 

kbar  in  the  wet  shot.  The  gage  in  C-7  on  the  back  surfaces  gave  a  peak 

stress  of  63  kbar  in  the  dry  shot  and  82  kbar  for  the  wet  shot.  An 
impedance  match  solution  using  the  previously  determined  Hugoniot 
equations  of  state  for  Coconino  sandstone  (Gregson  et  al.,  1963)  and 
for  C-7  (Keough,  1968)  gave  approximate  peak  stresses  of  7.1  kbar  at  the 
back  surface  for  the  dry  shot  and  106  kbar  for  the  wet  shot.  Waveforms 
at  both  the  internal  and  back  surface  gages  did  not  have  flat  tops  as 
expected  from  plate  impact  but  showed  stress  decay  immediately  behind  an 
early  peak.  This  waveform  is  characteristic  of  a  shock  being  attenuated 
by  and  overtaken  by  rarefaction.  This  interpretation  checks  with  the 
apparent  attenuation  between  the  in-material  gages  and  the  back  surface 
gages,  but  the  rarefaction  from  the  back  of  the  1/2-inch  aluminum  flyer 
could  not  overtake  the  shot  at  1/2-inch  depth  in  the  dry  sandstone.  Our 
only  explanation  at  present  is  that  a  stress  relaxation  occurs  after  the 
plate  impact  on  the  dry  sandstone.  The  internal  gages  in  the  wet  sand¬ 
stone  did  not  record  long  enough  to  show  whether  significant  stress  decay 
occurred  after  the  initial  peaks.  The  C-7  gage  on  the  back  of  th 3  wet 
shot  showed  a  decaying  stress. 


The  next  shots  fired  in  Coconino  sandstone,  one  dry  and  one  water 
saturated,  were  similarly  impacted  with  a  1/2-inch  Al  flyer  plate  driven 
by  a  P-80  plane  wave  lens  and  2  inches  of  Baratol.  Each  shot  was  instru¬ 
mented  with  manganic  stress  gages  at  the  impact  surface,  and  at  nominal 
distances  of  1/2,  1,  2,  and  :i  inches  from  the  impact  surface.  Each  gage 
was  cemented  between  two  sheets  of  Epon  0.003  inch  thick.  This  gage- 
sandwich  (total  thickness  about  0.007  inch)  was  then  placed  between  slabs 
of  sandstone.  In  the  case  of  the  water  saturated  shot,  the  slab  faces 
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were  coated  with  Homolite  epoxy.  The  penetration  of  the  epoxy  into  the 
sandstone  was  about  0.015  inch  or  about  2  grain  diameters.  The  sandstone 
slabs  were  then  saturated  by  allowing  water  to  enter  under  vacuum  from 
the  ends  of  each  slab.  Unfortunately,  in  the  water  saturated  shot,  the 
gage  leads  were  cut  by  the  explosive  chain,  and  measurements  were  not 
obtained.  In  the  dry  shot,  all  the  stress  gages  functioned  properly 
(see  Fig.  19)  except  for  the  one  on  the  impact  surface.  The  impact 
surface  gage  was  insulated  with  0.002-inch  of  Kapton  sheet,*  but  the 
gage  shorted  at  the  instant  of  impact.  Stress  amplitudes  and  rise  times 
for  the  four  successful  gages  are  given  in  Table  V.  Average  shock  veloc¬ 
ities  between  positions  are  also  given. 

The  rarefaction  from  the  back  side  of  the  flyer  plate  should  catch 
the  shock  front  at  about  the  position  of  the  third  gage.  The  exact  time 
of  catch-up  is  uncertain  because  the  rarefaction  velocity  in  the  com¬ 
pressed  sandstone  is  not  known.  To  estimate  the  time  of  catch-up,  we 
assumed  7.0  *ra/Hsec,  which  is  about  twice  the  shock  velocity.  Our 
previous  work  on  porous  playa  alluvium  indicates  that  such  a  value  is 
reasonable  (Petersen  et  al.,  1970). 

I I I  RECOMMENDATIONS 

The  four  sandstone  experiments  were  disappointing  in  that  the  few 
results  obtained  cannot  be  given  a  clear-cut  interpretation.  At  some 
future  date,  a  program  to  study  stress  wave  propagation  through  thick 
specimens  of  porous  rock  would  be  desirable.  This  would  necessarily 
be  a  larger  effort  with  additional  work  in  gage  development  and  emplace¬ 
ment  techniques  to  ensure  useful  results.  With  further  development  for 
use  in  explosive  experiments,  particle  velocity  gages  (Petersen  et  al. , 

1970)  may  prove  to  bi  superior  to  stress  gages  for  use  in  highly 
porous  rocks . 


LTenCe  Laboratory 

several  hundred  kblr  insulate  manganin  gages  in  aluminum  at 
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RECORO  FROM  GAGE  NO.  4 


FIGURE  19  GAGE  RECORDS  FOR  DRY  COO 
3  inches  from  impact. 


GP-7942-4 


RECORD  FROM  GAGE  NO.  5 


SANDSTONE.  Gages  at  1/2,  1,  2,  and 


Table  V 

DRY  COCONINO  SANDSTONE  GAGE  RECORDS 


Gage 

Gage  Distance 

Shock 

Stress 

Average  Shock 
Velocity 

No. 

from  Impact 

Rise-Time 

Amplitude 

Comments 

Between  Gages 

(mm) 

(p,sec) 

(kbar) 

(mm/pisec) 

1 

0 

— 

— 

Gage  shorted 
at  impact 

3.69 

2 

13.76 

s(e 

0.43 

114 

Peak 

74 

Level  position 

30 

Minimum 

52 

Final  level 
position 

3.64 

3 

28.34 

0.19 

90 

Peak 

45 

Level  position 

2.69 

4 

53.92 

0.25 

31 

2.32 

5 

78.30 

0.55 

26 

Peak 

16 

Level  position 

*  Structure  seen  in  rise. 
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